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Abstract. Wc study space-time fluctuations around a characteristic line for a one- 
dimensional interacting system known as the random average process. The state 
of this system is a real- valued function on the integers. New values of the function 
are created by averaging previous values with random weights. The fluctuations 
analyzed occur on the scale n^^^ where n is the ratio of macroscopic and microscopic 
scales in the system. The limits of the fluctuations are described by a family of 
Gaussian processes. In cases of known product-form invariant distributions, this 
limit is a two-parameter process whose time marginals are fractional Brownian 
motions with Hurst parameter j. Along the way we study the limits of quenched 
mean processes for a random walk in a space-time random environment. These 
limits also happen at scale n^^'^ and are described by certain Gaussian processes 
that we identify. In particular, when we look at a backward quenched mean process, 
the limit process is the solution of a stochastic heat equation. 



1. Introduction 

Fluctuations for asymmetric interacting systems. An asymmetric interacting 
system is a random process = {^"^-(/c) : /c G /C} of many components a"r(^) tliat 
influence each otliers' evolution. Asymmetry means here that the components have an 
average drift in some spatial direction. Such processes are called interacting particle 
systems because often these components can be thought of as particles. 

To orient the reader, let us first think of a single random walk {Xr : r = 0, 1, 2, . . . } 
that evolves by itself. For random walk we scale both space and time by n because on 
this scale we see the long-term velocity: n'^X^nt] ^ tv as n ^ oo where v = EXi. 
The random walk is diffusive which means that its fluctuations occur on the scale 
n^/^, as revealed by the classical central limit theorem: n~^^'^(X^nt} — ntv) converges 
weakly to a Gaussian distribution. The Gaussian limit is universal here because it 
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arises regardless of the choice of step distribution for the random walk, as long as a 
square-integrability hypothesis is satisfied. 

For asymmetric interacting systems we typically also scale time and space by the 
same factor n, and this is known as Euler scaling. However, in certain classes of one- 
dimensional asymmetric interacting systems the random evolution produces fluctua- 
tions of smaller order than the natural diffusive scale. Two types of such phenomena 
have been discovered. 

(i) In Hammersley's process, in asymmetric exclusion, and in some other closely 
related systems, dynamical fluctuations occur on the scale n^^^. Currently known 
rigorous results suggest that the Tracy- Widom distributions from random matrix 
theory are the universal limits of these n^^^ fluctuations. 

The seminal works in this context are by Baik, Deift and Johansson |2] on Ham- 
mersley's process and by Johansson on the exclusion process. We should point 
out though that j3] does not explicitly discuss Hammersley's process, but instead the 
maximal number of planar Poisson points on an increasing path in a rectangle. One 
can intrepret the results in [3] as fluctuation results for Hammersley's process with a 
special initial configuration. The connection between the increasing path model and 
Hammersley's process goes back to Hammersley's paper ^H]- It was first utilized by 
Aldous and Diaconis (who also named the process), and then further in the papers 

[23 123- 

(ii) The second type has fluctuations of the order n^^^ and limits described by a 
family of self-similar Gaussian processes that includes fractional Brownian motion 
with Hurst parameter |. This result was first proved for a system of independent 
random walks [201 • One of the main results of the current paper shows that the n^^^ 
fluctuations also appear in a family of interacting systems called random average pro- 
cesses in one dimension. The same family of limiting Gaussian processes appears here 
too, suggesting that these limits are universal for some class of interacting systems. 

The random average processes (RAP) studied in the present paper describe a ran- 
dom real-valued function on the integers whose values evolve by jumping to random 
convex combinations of values in a finite neighborhood. It could be thought of as a 
caricature model for an interface between two phases on the plane, hence we call the 
state a height function. RAP is related to the so-called linear systems discussed in 
Chapter IX of Liggett's monograph [221 • RAP was introduced by Ferrari and Pontes 
P!lj who studied the fluctuations from initial linear slopes. In particular, they discov- 
ered that the height over the origin satisfies a central limit theorem in the time scale 
t^/^. The Ferrari-Fontes results suggested RAP to us as a fruitful place to investigate 
whether the n^/^ fluctuation picture discovered in jSHI for independent walks had any 
claim to universality. 

There are two ways to see the lower order dynamical fluctuations. 
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(1) One can take deterministic initial conditions so that only dynamical random- 
ness is present. 

(2) Even if the initial state is random with central limit scale fluctuations, one 
can flnd the lower order fluctuations by looking at the evolution of the process 
along a characteristic curve. 

Articles j2] and studied the evolutions of special deterministic initial states of 
Hammersley's process and the exclusion process. Recently Ferrari and Spohn ^3] 
have extended this analysis to the fluctuations across a characteristic in a stationary 
exclusion process. The general nonequilibrium hydrodynamic limit situation is still 
out of reach for these models. [20] contains a tail bound for Hammersley's process 
that suggests n^/^ scaling also in the nonequilibrium situation, including along a 
shock which can be regarded as a "generalized" characteristic. 

Our results for the random average process are for the general hydrodynamic limit 
setting. The initial increments of the random height function are assumed indepen- 
dent and subject to some moment bounds. Their means and variances must vary 
sufficiently regularly to satisfy a Holder condition. Deterministic initial increments 
qualify here ClS db special case of independent. 

The classification of the systems mentioned above (Hammersley, exclusion, in- 
dependent walks, RAP) into n^^^ and n^^^ fiuctuations coincides with their clas- 
sification according to type of macroscopic equation. Independent particles and 
RAP are macroscopically governed by linear first-order partial differential equations 
Ut + hux = 0. In contrast, macroscopic evolutions of Hammersley's process and the 
exclusion process obey genuinely nonlinear Hamilton- Jacobi equations Ut + f{ux) = 
that create shocks. 

Suppose we start off one of these systems so that the initial state fluctuates on the 
n^/^ spatial scale, for example in an stationary distribution. Then the fluctuations of 
the entire system on the n^/^ scale simply consist of initial fluctuations transported 
along the deterministic characteristics of the macroscopic equation. This is a conse- 
quence of the lower order of dynamical fluctuations. When the macroscopic equation 
is linear this is the whole picture of diffusive fluctuations. In the nonlinear case the 
behavior at the shocks (where characteristics merge) also needs to be resolved. This 
has been done for the exclusion process [23] and for Hammersley's process |29j . 

Random walk in a space-time random environment. Analysis of the random 
average process utilizes a dual description in terms of backward random walks in a 
space-time random environment. Investigation of the fluctuations of RAP leads to a 
study of fluctuations of these random walks, both quenched invariance principles for 
the walk itself and limits for the quenched mean process. The quenched invariance 
principles have been reported elsewhere |21]. The results for the quenched mean 
process are included in the present paper because they are intimately connected to 
the random average process results. 
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We look at two types of processes of quenched means. We call them forward and 
backward. In the forward case the initial point of the walk is fixed, and the walk 
runs for a specified amount of time on the space-time lattice. In the backward case 
the initial point moves along a characteristic, and the walk runs until it reaches the 
horizontal axis. Furthermore, in both cases we let the starting point vary horizontally 
(spatially), and so we have a space-time process. In both cases we describe a limiting 
Gaussian process, when space is scaled by n^/^, time by n, and the magnitude of the 
fluctuations by n^/^. In particular, in the backward case we find a limit process that 
solves the stochastic heat equation. 

There are two earlier papers on the quenched mean of this random walk in a space- 
time random environment. These previous results were proved under assumptions 
of small enough noise and finitely many possible values for the random probabili- 
ties. Bernabei jH] showed that the centered quenched mean, normalized by its own 
standard deviation, converges to a normal variable. Then separately he showed that 
this standard deviation is bounded above and below on the order n^/^. Bernabei 
has results also in dimension 2, and also for the quenched covariance of the walk. 
Boldrighini and Pellegrinotti [H] also proved a normal limit in the scale n^^^ for what 
they term the "correction" caused by the random environment on the mean of a test 
function. 

Finite-dimensional versus process-level convergence. Our main results all 
state that the finite-dimensional distributions of a process of interest converge to 
the finite-dimensional distributions of a certain Gaussian process specified by its co- 
variance function. We have not proved process-level tightness, except in the case of 
forward quenched means for the random walks where we compute a bound on the 
sixth moment of the process increment. 

Further relevant literature. It is not clear what exactly are the systems "closely 
related" to Hammersley's process or exclusion process, alluded to in the beginning 
of the Introduction, that share the n^^^ fluctuations and Tracy- Widom limits. The 
processes for which rigorous proofs exist all have an underlying representation in 
terms of a last-passage percolation model. Another such example is "oriented digital 
boiling" studied by Gravner, Tracy and Widom JH] • (This model was studied earlier 
in j27j and [201 under different names.) 

Fluctuations of the current were initially studied from the perspective of a moving 
observer traveling with a general speed. The fluctuations are diffusive, and the lim- 
iting variance is a function of the speed of the observer. The special nature of the 
characteristic speed manifests itself in the vanishing of the limiting variance on this 
diffusive scale. The early paper of Ferrari and Fontes |T2] treated the asymmetric 
exclusion process. Their work was extended by Balazs j3] to a class of deposition 
models that includes the much-studied zero range process and a generalization called 
the bricklayers' process. 
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Work on the fluctuations of Hammersley's process and tlie exclusion process has 
connections to several parts of mathematics. Overviews of some of these links appear 
in papers j2], JHI; and J7]. General treatments of large scale behavior of interacting 
random systems can be found in j^I], [22], [22], [22], and [22] • 

Organization of the paper. We begin with the description of the random average 
process and the limit theorem for it in Section |2j Section |2] describes the random 
walk in a space-time random environment and the limit theorems for quenched mean 
processes. The proofs begin with Section E) that lays out some preliminary facts on 
random walks. Sections and |B] prove the fluctuation results for random walk, and 
the final Section [3 proves the limit theorem for RAP. 

The reader only interested in the random walk can read Section |2] and the proofs 
for the random walk limits independently of the rest of the paper, except for certain 
definitions and a hypothesis which have been labeled. The RAP results can be read 
independently of the random walk, but their proofs depend on the random walk 
results. 

Notation. We summarize here some notation and conventions for quick reference. 
The set of natural numbers is N = {1, 2, 3, . . . }, while Z.,. = {0, 1, 2, 3, ... } and 
= [0, cxd). On the two dimensional integer lattice l? standard basis vectors are 
ei = (1,0) and 62 = (0, 1). The e2-direction represents time. 

We need several different probability measures and corresponding expectation op- 
erators. P (with expectation E) is the probability measure on the space f2 of en- 
vironments 10. P is an i.i.d. product measure across the coordinates indexed by the 
space-time lattice 1? . P (with expectation E) is the probability measure of the initial 
state of the random average process. E"^ is used to emphasize that an expectation 
over initial states is taken with a fixed environment lo. Jointly the environment and 
initial state are independent, so the joint measure is the product P(g)P. (with ex- 
pectation E^^ is the quenched path measure of the random walks in environment lo. 
The annealed measure for the walks is P = / P^^idiS). Additionally, we use P and 
E for generic probability measures and expectations for processes that are not part 
of this specific set-up, such as Brownian motions and limiting Gaussian processes. 

The environments 0; G f2 are configurations lo = {ux^t '■ [x, r) G 7?) of vectors 
indexed by the space-time lattice 7?. Each element uJx^t is a probability vector of 
length 2M + 1, denoted also by Ut-{x) = uJx,t, and in terms of coordinates u^(x) = 
{ur{x,y) : — M < y < M). The environment at a fixed time value r is o;^ = {uj^^t '■ 
X E Z). Translations on fl are defined by {Tx^T^^)y^s = ^x+y,T+s- 

[x\ = max{?7, G Z : 77, < x} is the lower integer part of a real x. Throughout, 
C denotes a constant whose exact value is immaterial and can change from line 
to line. The density and cumulative distribution function of the centered Gaussian 
distribution with variance o"^ are denoted by if„2[x) and ^^2[x). {B{t) : t > 0} 
is one-dimensional standard Brownian motion, in other words the Gaussian process 
with covariance EB{s)B{t) = s At. 



6 



MARTON BALAZS, FIRAS RASSOUL-AGHA, AND TIMO SEPPALAINEN 



Acknowledgements. The authors thank P. Ferrari and L. R. Pontes for comments 
on article [E] and J. Swanson for helpful discussions. 

2. The random average process 

The state of the random average process (RAP) is a height function o" : Z — >■ M. 
It can also be thought of as a sequence o = {cr{i) : i E Z) E M.^ where a{i) is the 
height of an interface above site i. The state evolves in discrete time according to the 
following rule. At each time point r = 1, 2, 3, . . . and at each site /c G Z, a random 
probability vector Ur{k) = {ur{k,j) : —M < j < M) of length 2M + 1 is drawn. 
Given the state (Jt-_i = {(Tr^i^i) : i G Z) at time r — 1, the height value at site k is 
then updated to 

(^Ak)= ^ Ur{k, j)ar^i{k + j). (2.1) 

j--\j\<M 

This update is performed independently at each site k to form the state o"r = (cr^(A;) : 
A; G Z) at time r. The same step is repeated at the next time r + 1 with new 
independent draws of the probability vectors. 

So, given an initial state (Tq, the process ar is constructed with a collection {u^{k) : 
r G N, A; G Z} of independent and identically distributed random vectors. These 
random vectors are defined on a probability space (fi, ©,P). If (Tq is also random 
with distribution P, then (Jq and the vectors {^^-(A;)} are independent, in other words 
the joint distribution is P ® P. We write u^{k) to make explicit the dependence on 
Co" G fi. E will denote expectation under the measure P. M is the range and is a fixed 
finite parameter of the model. P-almost surely each random vector Ur{k) satisfies 

M 

< Ur{k,j) < 1 for all -M <j<M, and ^ u^ikj) = 1. 

j=-M 

It is often convenient to allow values M^(fc, j) for all j. Then automatically ^^-(fc, j) = 
for |j| > M. 
Let 

p(0,j)=Emo(0,j) 

denote the averaged probabilities. Throughout the paper we make two fundamental 
assumptions. 

(i) Pirst, there is no integer h > 1 such that, for some x G Z, 

^p(0,x + kh) = 1. 

fcez 

This is also expressed by saying that the span of the random walk with jump 
probabilities p{0,j) is 1 [TTJ page 129]. It follows that the group generated by 
{x G Z : p(0, x) > 0} is all of Z, in other words this walk is aperiodic in Spitzer's 
terminology ^T] . 



RANDOM AVERAGE PROCESS 7 

(ii) Second, we assume that 

P{maxMo(0,j) < 1} > 0. (2.2) 

i 

If this assumption fails, then P-almost surely for each (/c, r) there exists j = j{k,T) 
such that Ur{k,j) = 1. No averaging happens, but instead ar{k) adopts the value 
ar-i{k + j). The behavior is then different from that described by our results. 

No further hypotheses are required of the distribution P on the probability vec- 
tors. Deterministic weights u'^{k,j) = p{0,j) are also admissible, in which case ()2.2p 
requires maxjp(0, j) < 1. 

In addition to the height process we also consider the increment process rjj. = 
{r]r{i) : i E Z,) defined by 

rjrii) = a-r{i) - a^ii - 1). 

From 1)2.11) one can deduce a similar linear equation for the evolution of the increment 
process. However, the weights are not necessarily nonnegative, and even if they are, 
they do not necessarily sum to one. 

Next we define several constants that appear in the results. 

D{uj) = ^xu'^{0,x) (2.3) 

is the drift at the origin. Its mean is V = K{D) and variance 

al = E[{D-Vn (2.4) 
A variance under averaged probabilities is computed by 

al = ^ix-V)'piO,x). (2.5) 

Define random and averaged characteristic functions by 

(f)^{t) = ^u^(0,x)e**^ and (f)a{t) = E0'^(t) = J]p(0,x)e^*^ (2.6) 

x&Z x&L 

and then further 

A(t)=E[|r(t)P] and A(t) = |0,(t)|^ (2.7) 
Finally, define a positive constant /5 by 

^^ri-A^ (2.8) 

The assumption of span 1 implies that |0a(^)| = 1 only at multiples of 27r. Hence the 
integrand above is positive at t 7^ 0. Separately one can check that the integrand has 
a finite limit as t — >■ 0. Thus /9 is well-defined and finite. 

In Section m we can give these constants, especially /3, more probabilistic meaning 
from the perspective of the underlying random walk in random environment. 
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For the limit theorems we consider a sequence a" of the random average processes, 
indexed by n G N = {1, 2, 3, . . . }. Initially we set ctq (0) = 0. For each n we assume 
that the initial increments {f]Q{i) : i & 1j} are independent random variables, with 

E[r]^{i)] = Q{i/n) and Var[r/^(z)] = v{i/n). (2.9) 

The functions g and v that appear above are assumed to be uniformly bounded 
functions on M and to satisfy this local Holder continuity: 

For each compact interval [a, 6] C M there exist 

C = C{a, b) < oo and 7 = 7(0, b) > 1/2 such that (2.10) 

\g{x) — g{y)\ + \v{x) — v{y)\ < C \x — y\'^ for x,y E [a, b]. 

The function v must be nonnegative, but the sign of g is not restricted. Both functions 
are allowed to vanish. In particular, our hypotheses permit deterministic initial 
heights which implies that v vanishes identically. 

The distribution on initial heights and increments described above is denoted by 
P. We make this uniform moment hypothesis on the increments: 

there exists a > such that sup E[|r/^(i)|2+"] < 00. (2.11) 

We assume that the processes a" are all defined on the same probability space. The 
environments u that drive the dynamics are independent of the initial states {(Tq}, so 
the joint distribution of (cj, {ctq}) is P ® P. When computing an expectation under 
a fixed u we write E*^. 

On the larger space and time scale the height function is simply rigidly translated 
at speed b = —V, and the same is also true of the central limit fluctuations of the 
initial height function. Precisely speaking, define a function U on M by U{0) = and 
U'{x) = g{x). Let {x,t) G M x ]R+. The assumptions made thus far imply that both 

n-Vf„,j(M)^f/(x-6t) (2.12) 

and 

^"ntj(M) -^^(3^-^^) ag([nxj - [nbt\)-nU{x-bt) 

1= 1= ^ U [2.16) 

in probability, as n 00. (We will not give a proof. This follows from easier 
versions of the estimates in the paper.) Limit ()2.12j) is the "hydrodynamic limit" of 
the process. The large scale evolution of the height process is thus governed by the 
linear transport equation 

Wt + bwx = 0. 

This equation is uniquely solved by w{x,t) = U{x — bt) given the initial function 
w(x, 0) = U{x). The lines x{t) = x + bt are the characteristics of this equation, 
the curves along which the equation carries information. Limit ()2.13p says that 
fluctuations on the diffusive scale do not include any randomness from the evolution, 
only a translation of initial fluctuations along characteristics. 
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We find interesting height fluctuations along a macroscopic characteristic hne 
x{t) = y + bt, and around such a hne on the microscopic spatial scale y/n. The 
magnitude of these fluctuations is of the order n^^^, so we study the process 

Zr.{t,r) = n-'/'{af^,^{[ny\ + [rV^\ + [ntb\) - a^i[ny\ + [rv^J)}, 

indexed by {t, r) G M+ x M, for a flxed ?/ G M. In terms of the increment process 77^, 
Zn{t, 0) is the net flow from right to left across the discrete characteristic lny\ + [risfej , 
during the time interval < s <t. 

Next we describe the limit of Zn- Recall the constants deflned in (|2.4j) . (|2.5p . and 
Combine them into a new constant 



~ o 2- (2-14) 

Let {B{t) : t > 0} be one-dimensional standard Brownian motion. Deflne two 
functions Tq and Tq on (M+ x M) x (M+ x R): 

r,((., g), (t, r)) = ^ r^*^'^ exp|-^(g - rf] dv (2.15) 



and 

ro((s,g),(t,r)) 



(2.16) 



/■oo 

/ P[aaB{s) > X- q]P[aaB{t) > x-r]dx 

JqVr 

- {l{r>,} I P[aaB{s) >x- q]P[aaB{t) <x-r]dx 

+ 1|,>,} J' P[aaB{s) <x- q]P[aaB{t) > x - r] rfx} 

/gAr 
^ka5(s) < X - q]P[aaB{t) <x-r] dx. 
■00 

The boundary values are such that rq((s, g), (t, r)) = ro((s, g), (t, r)) = if either 
s = or t = 0. We will see later that Tq is the limiting covariance of the backward 
quenched mean process of a related random walk in random environment. Fq is 
the covariance for fluctuations contributed by the initial increments of the random 
average process. (Hence the subscripts q for quenched and for initial time. The 
subscript on Vq has nothing to do with the argument (s, q).) 

The integral expressions above are the form in which Tq and Tq appear in the 
proofs. For Tq the key point is the limit ()5.19p which is evaluated earlier in ()4.5p . Fq 
arises in Proposition 17. 11 

Here are alternative succinct representations for Tq and Fq. Denote the centered 
Gaussian density with variance cr^ by 

^..(x) = -^exp{-^x^} (2.17) 
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and its distribution function by (x) = Lp„2 (y) dy. Then define 

^^2(x) = a'^ifa^ix) - x{l - $^2(x)), 
which is an antiderivative of $0-2 (x) — 1. In these terms, 

Tg{{s, q), (t, r)) = K^^2(t+s){\q -r\)- K^'^2|t_,| (|g - r|). 

and 

ro((s, g), (t, r)) = - r|) + ^.2^ (|g - r|) - ^.2(i+,) (|g - r|). 

Theorem 2.1. Assume (j2.2p anc? t/iat t/ie averaged probabilities p{0,j) = EMQ(0,j) 
/iai'e lattice span 1. Let g and v be two uniformly hounded functions on R that satisfy 
the local Holder condition ()2.1()p . For each n, let be a random average process 
normalized by cTq (0) = and whose initial increments {rj^li) : i G Z} are independent 
and satisfy ()2.9|) and ()2.11|) . Assume the environments uj independent of the initial 
heights {ctq : G N}. 

Fix y e M. Under the above assumptions the finite- dimensional distributions of the 
process {z„(t,r) : (t,r) G M+xM} converge weakly as n ^ 00 to the finite- dimensional 
distributions of the mean zero Gaussian process {z{t, r) : (t, r) G M_|_ x M} specified 
by the covariance 

Ez{s, q)z{t, r) = f?(y)'r,((s, g), (t, r)) + v{y)To{{s, q), {t, r)). (2.18) 



The statement means that, given space-time points (ti, ri), . . . , (t^, r^), the M*^- 
valued random vector (-2,1(^1, ri), . . . , Zn(tk, r^)) converges in distribution to the ran- 
dom vector ri), . . . , 2;(tfc, rfc)) as — > 00. The theorem is also vahd in cases 
where one source of randomness has been turned off: if initial increments around 
\ny\ are deterministic then v{y) = 0, while if D{uj) = V then o"^ = 0. The case 
cr|, = contains as special case the one with deterministic weights u^{k,j) = p{0,j). 

If we consider only temporal correlations with a fixed r, the formula for the covari- 
ance is as follows: 

Ez{s, r)z{t, r) = ^g{yY{y/JTt - y/t^S ) 

a (2-19) 
+ -^v{y) (v^ + - + t ) fors<t. 



/27r 

Remark 2.1. The covariances are central to our proofs but they do not illuminate 
the behavior of the process z. Here is a stochastic integral representation of the 
Gaussian process with covariance ()2.18|1 : 



z(t, r) = g{y)(Tay^ / / (/)^2(j„,)(r - x) dW{s, x) 

JJ[o,t]xm (2 20) 

+ V 'v{y) I sign(x — r)$o-2^ (— |x — r| ) (ii?(x). 
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Above ly is a two-parameter Brownian motion defined on ]R_|_ x M, 5 is a one- 
parameter Brownian motion defined on M, and W and B are independent of eacli 
otlier. Tlie first integral represents the space-time noise created by the dynamics, 
and the second integral represents the initial noise propagated by the evolution. The 
equality in ()2.2m) is equality in distribution of processes. It can be verified by checking 
that the Gaussian process defined by the sum of the integrals has the covariance 

(inHi). 

One can readily see the second integral in (|2.2Up arise as a sum in the proof. It is 
the limit of F"(t,r) defined below equation (j7.1|) . 

One can also check that the right-hand side of (j2.2m) is a weak solution of a 
stochastic heat equation with two independent sources of noise: 

zt = \<tI Zrr + Q{:y)<ya^/^W + l^^al B", z{0, r) = 0. (2.21) 

W is space-time white noise generated by the dynamics and B" the second derivative 
of the one- dimensional Brownian motion that represents initial noise. This equation 
has to be interpreted in a weak sense through integration against smooth compactly 
supported test functions. We make a related remark below in Section 13.21 for limit 
processes of quenched means of space-time RWRE. 

The simplest RAP dynamics averages only two neighboring height values. By 
translating the indices, we can assume that p(0, —1) + p(0, 0) = 1. In this case the 
evolution of increments is given by the equation 

r]^{k) = Ur{k, 0)r]r-i{k) + n,(A; - 1, -l)r/,_i(A; - 1). (2.22) 

There is a queueing interpretation of sorts for this evolution. Suppose rjr-iik) denotes 
the amount of work that remains at station k at the end of cycle r — 1. Then during 
cycle r, the fraction m,-(/c, —1) of this work is completed and moves on to station 
k + 1, while the remaining fraction UT-{k^ 0) stays at station k for further processing. 

In this case we can explicitly evaluate the constant (3 in terms of the other quanti- 
ties. In a particular stationary situation we can also identify the temporal marginal 
of z in (j2.19p as a familiar process. (A probability distribution /i on the space is 
an invariant distribution for the increment process if it is the case that when rjQ has 
/i distribution, so does r]r for all times r G 2+.) 

Proposition 2.2. Assume p(0, -1) + p(0, 0) = 1. 

(a) Then 

/?= ^E[mo(0,0)mo(0,-1)]. (2.23) 

(b) Suppose further that the increment process rj^- possesses an invariant distribu- 
tion fi in which the variables {r]{i) : i E Z} are i.i.d. with common mean g = E'^[ri{i)] 
and variance v = E^[ri[i)'^] — . Then v = ng^ . 
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Suppose that in Theorem \2.1\ each rj'^ = rj^- is a stationary process with marginal /i. 
Then the limit process z has covariance 

Ez{s, q)z{t, r) = kq" (^.g, (|g - r|) + - r |) - *.2|,„,| (|g - r |)) . (2.24) 

In particular, for a fixed r the process {z{t,r) : t G M+} has covariance 

2 

Ezis,r)z{t,r) = ^^{^s + Vt - ^/\t^\) . (2.25) 
V 27r 

In other words, process z{-,r) is fractional Brownian motion with Hurst parameter 
1/4. 

To rephrase the connection ()2.24|1 - ()2.25|) . the process {z(t,r)} in ()2.24|1 is a cer- 
tain two-parameter process whose marginals along the first parameter direction are 
fractional Brownian motions. 

Ferrari and Pontes J3] showed that given any slope p, the process r]r started from 
deterministic increments rio{x) = px converges weakly to an invariant distribution. 
But as is typical for interacting systems, there is little information about the invariant 
distributions in the general case. The next example gives a family of processes and 
i.i.d. invariant distributions to show that part (b) of Proposition 12.21 is not vacuous. 
Presently we are not aware of other explictly known invariant distributions for RAP. 

Example 2.1. Fix integer parameters m > j > 0. Let {Ur{k, —1) : r E N, k E Z} he 
i.i.d. beta-distributed random variables with density 

(j-l)!(m-j-l)!" 

on (0, 1). Set Ur{k, 0) = 1 — u^-^k, —1). Consider the evolution defined by ()2.22p with 
these weights. Then a family of invariant distributions for the increment process 
rj^ = {rir{k) : k E Tj) is obtained by letting the variables {ri{k)} be i.i.d. gamma 
distributed with common density 

fix) = -^—-\e-'^i\xr-' (2.26) 
[m — 1)1 

on M+. The family of invariant distributions is parametrized by < A < oo. Under 
this distribution E[?7(A;)] = m/X and Var[?7(/c)] = m/A^. 

One motivation for the present work was to investigate whether the limits found 
in [201 for fluctuations along a characteristic for independent walks are instances of 
some universal behavior. The present results are in agreement with those obtained 
for independent walks. The common scaling is n^^'^. In that paper only the case 
r = of Theorem 12.11 was studied. For both independent walks and RAP the limit 
z{- ,0) is a mean- zero Gaussian process with covariance of the type 

Ez{s, 0)z{t, 0) = Ci {Vs + t - y/t-S ) + C2 (a/s + Vs + t ) 
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where ci is determined by the mean increment and C2 by the variance of the increment 
locally around the initial point of the characteristic. Furthermore, as in Proposition 
I2.2f b). for independent walks the limit process specializes to fractional Brownian 
motion if the increment process is stationary. 

These and other related results suggest several avenues of inquiry. In the introduc- 
tion we contrasted this picture of "n}!^ fluctuations and fractional Brownian motion 
limits with the n^/^ fluctuations and Tracy- Widom limits found in exclusion and 
Hammersley processes. Obviously more classes of processes should be investigated 
to understand better the demarcation between these two types. Also, there might be 
further classes with different limits. 

Above we assumed independent increments at time zero. It would be of interest to 
see if relaxing this assumption leads to a change in the second part of the covariance 
fl2.18|) . [The first part comes from the random walks in the dual description and 
would not be affected by the initial conditions.] However, without knowledge of some 
explicit invariant distributions it is not clear what types of initial increment processes 
{?7o(/c)} are worth considering. Unfortunately finding explicit invariant distributions 
for interacting systems seems often a matter of good fortune. 

We conclude this section with the dual description of RAP which leads us to study 
random walks in a space-time random environment. Given a;, let {X.^ '■ s G 
denote a random walk on Z that starts at Xq^ = i, and whose transition probabilities 
are given by 

P^iXl';, = y I Xy = x)= <_,(a:, y - x). (2.27) 

P'^ is the path measure of the walk X*'"^, with expectation denoted by E'^. Compar- 
ison of and gives 

= ^P-(Xr = J I X^'^ = ^)a,.^U) = . (2.28) 

j 

Iteration and the Markov property of the walks X*'"^ then lead to 

ar{i)=E''[ao{Xl'-)]. (2.29) 

Note that the initial height function o"o is a constant under the expectation E'^. 

Let us add another coordinate to keep track of time and write X*'"^ = (X*''^, r — s) 
for s > 0. Then X*'"^ is a random walk on the planar lattice that always moves 
down one step in the e2-direction, and if its current position is (x, n), the ei-coordinate 
of its next position is x + y with probability Un{x,y). We shall call it the backward 
random walk in a random environment. In the next section we discuss this walk and 
its forward counterpart. 

3. Random walk in a space-time random environment 

3.1. Definition of the model. We consider here a particular random walk in ran- 
dom environment (RWRE). The walk evolves on the planar integer lattice Z^, which 
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we think of as space-time: the first component represents one-dimensional discrete 
space, and the second represents discrete time. We denote by 62 the unit vector in 
the time-direction. The walks will not be random in the e2-direction, but only in the 
spatial Ci-direction. 

We consider forward walks and backward walks X^^. The subscript m G Z+ 
is the time parameter of the walk and superscripts are initial points: 

Z^'- = X^'- = (z,r)GZ2. (3.1) 

The forward walks move deterministically up in time, while the backward walks move 
deterministically down in time: 

Z^- = (Z^-, T + m) and X^" = (X^", r - m) for m > 0. 

Since the time components of the walks are deterministic, only the spatial components 
Zli[ and are really relevant. We impose a finite range on the steps of the walks: 
there is a fixed constant M such that 

|zr+i-^;;1<M and \xi;,l,-Xl;:\<M. (3.2) 

A note of advance justification for the setting: The backward walks are the ones 
relevant to the random average process. Distributions of forward and backward walks 
are obvious mappings of each other. However, we will be interested in the quenched 
mean processes of the walks as we vary the final time for the forward walk or the 
initial space-time point for the backward walk. The results for the forward walk form 
an interesting point of comparison to the backward walk, even though they will not 
be used to analyze the random average process. 

An environment is a configuration of probability vectors uj = (ur{x) : (x, r) G 
where each vector Ur{x) = {ur{x,y) : —M < y < M) satisfies 

M 

< Ur{x, y) < 1 for all -M < y < M, and ^ m^(x, y) = 1. 

An environment u is a. sample point of the probability space {Q, (5,P). The sample 
space is the product space Q = V^'^ where V is the space of probability vectors 
of length 2M + 1, and & is the product cx-field on Q induced by the Borel sets 
on V. Throughout, we assume that P is a product probability measure on Q such 
that the vectors {ur{x) : {x,t) G Z^} are independent and identically distributed. 
Expectation under P is denoted by E. When for notational convenience we wish to 
think of Ur{x) as an infinite vector, then Ur{x,y) = for \y\ > M. We write u^{x,y) 
to make explicit the environment u, and also uJx,t = Ur{x) for the environment at 
space-time point (x, r). 

Fix an environment uj and an initial point (i, r). The forward and backward walks 
Zl^ and X'l^ (m > 0) are defined as canonical Z^- valued Markov chains on their 
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path spaces under the measure P'^ determined by the conditions 

P-{Z^'^ = i^, r)} = l, 

P-{Z:V^l = {y,T + S + l)\Zl'- = {X,T + S)}= Ur+s{x, y - x) 

for the forward walk, and by 

P-{Xr = r)] = 1, 
P^^^+i = {y,r-s-l)\Xl'^ = {x,r-s)} = y - x) 

for the backward walk. By dropping the time components r, r ± s and r ± s ± 1 
from the equations we get the corresponding properties for the spatial walks Z^''^ 
and Xl''^. When we consider many walks under a common environment u, it will 
be notationally convenient to attach the initial point {i, r) to the walk and only the 
environment u to the measure P'^. 

P'^ is called the quenched distribution, and expectation under P'^ is denoted by E^. 
The annealed distribution and expectation are P(-) = EP'^(-) and E{-) = KE'^{-). 
Under P both and are ordinary homogeneous random walks on Z with jump 
probabilities p{i, i + j) = p{0, j) = Emo(0,j). These walks satisfy the law of large 
numbers with velocity 

V = J2piO,j)j. (3.3) 

jez 

As for RAP, we also use the notation b = —V. 

3.2. Limits for quenched mean processes. We start by stating the quenched 
invariance principle for the space-time RWRE. {B{t) : t > 0} denotes standard one- 
dimensional Brownian motion. -Dr[0, oo) is the space of real- valued cadlag functions 
on [0, oo) with the standard Skorohod metric Recall the definition ()2.5|) of the 
variance of the annealed walk, and assumption ()2.2j) that guarantees that the 
quenched walk has stochastic noise. 

Theorem 3.1. [21] Assume ()2.2|) . Then for almost every u, under P^ the process 
j7,-i/2^j5^o^o^ — nt\^) converges weakly to the process B{a'^t) on the path space D^[0, oo) 
as n oo. 

Assume further that aj^ > so that the environment is not degenerate. Then 
we have these bounds on the variance of the quenched mean: there exist constants 
< Ci, 6*2 < oo such that for all n 

Cm^/^ < E[(E"(X°'°) - nVy] < Csn^/'. (3.4) 

Quite obviously, X°'° and Z^'^ are interchangeable in the above theorem. Bounds 
fl3.4p suggest the possibility of a weak limit for the quenched mean on the scale n^^^. 
Such results are the main point of this section. 
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For t > 0, r G M we define scaled, centered quenched mean processes 

a„(t,r) = n-''\E-{z[^f^'') - [r^\ - [nt\V] (3.5) 
for tlie forward walks, and 

for the backward walks. In words, the process an follows forward walks from level 
to level \nt\ and records centered quenched means. Process ?/.„ follows backward 
walks from level \nt\ down to level and records the centered quenched mean of the 
point it hits at level 0. The initial points of the backward walks are translated by 
the negative of the mean drift \nth\. This way the temporal processes a„(-,r) and 
yn{-,r) obtained by fixing r are meaningful processes. 
Random variable yn{t,r) is not exactly centered, for 

E|/„(t,r) = n-^^^{[ntb\ - [nt\b). (3.7) 

Of course this makes no difference to the limit. 

Next we describe the Gaussian limiting processes. Recall the constant k defined 
in and the function Tg defined in fITTHll . Let {a{t,r) : {t,r) e M+ x M} and 

{y{t,r) : {t,r) G IR+ x M} be the mean zero Gaussian processes with covariances 

Ea{s,q)a{t,r) = Tq{{s At,q),{s At,r)) 

and 

Ey{s, q)y{t, r) = r^((s, g), (t, r)) 

for s,t > and g, r G M. When one argument is fixed, the random function r h-s- y(t, r) 
is denoted by y{t, ■) and 1 1— > y{t, r) by y{-, r). From the covariances follows that at a 
fixed time level t the spatial processes a{t, ■) and y{t, ■) are equal in distribution. 
We record basic properties of these processes. 

Lemma 3.1. The process {y{t,r)} has a version with continuous paths as functions 
of (t,r). Furthermore, it has the following Markovian structure in time. Given = 
to < ^1 < ■ ■ ■ < ^n, let {y(ti — •) : 1 < i < n} be independent random functions 
such that y{ti — ti-i, ■) has the distribution of y(ti — ■) for i = 1, . . . ,n. Define 
y*{ti,r) = y{ti, r) for r G and then inductively for i = 2, . . . ,n and r G M, 

y*{U,r) = / v?a2(t.-ti_i)(^^)y*(i»-i,r + u)rfu + y(ti-ti_i,r). (3.8) 
Jr 

Then the joint distribution of the random functions {y*{ti, ■) : 1 <i <n} is the same 
as that of {y{ti, ■) '■ 1 < i < n} from the original process. 

Sketch of proof . Consider (s, q) and {t, r) varying in a compact set. From the covari- 
ance comes the estimate 

E[{y{s,q)-y{t,r)y] < C {\s - t\'/' + \q - r\) (3.9) 
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from which, since the integrand is Gaussian, 

E[{y{s, q) - y{t, r)f'] < C {\s - t\"^ + \q - r\f < C \\ (s, q) - (t, r)f^ . (3.10) 

Kolmogorov's criterion imphes the existence of a continuous version. 

For the second statement use ()3.8|1 to express a hnear combination Yll=i ^iy*(ti^ ''^i) 
in the form 



n n „ 

i=i i=i -^^ 



x) \i{dx) 



where the signed measures Aj are hnear combinations of Gaussian distributions. Use 
this representation to compute the variance of the hnear combination on the left-hand 
side (it is mean zero Gaussian). Observe that this variance equals 

Lemma 3.2. The process {a{t, r)} has a version with continuous paths as functions of 
{t, r). Furthermore, it has independent increments in time. A more precise statement 
follows. Given = to < ti < ■ ■ ■ < tn, let {d{ti — ■) : 1 < i < n} be independent 
random functions such that a(ti — tj-i,-) has the distribution of a(ti — ti-i,-) for 
2 = 1, . . . , n. Define a*(ti,r) = a(ti,r) forr G M, and then inductively for i = 2, . . . ,n 
and r G M, 



a*{ti,r) = a*{ti_i,r) + / ip^2t^_^{u)a{ti - ti_i,r + u) du. (3.11) 

Jr 

Then the joint distribution of the random functions {a*(tj, ■) : 1 < i < n} is the same 
as that of {a(tj, ■) '. 1 < i < n} from the original process. 

The proof of the lemma above is similar to the previous one so we omit it. 

Remark 3.1. Processes y and a have representations in terms of stochastic integrals. 
As in Remark 12. II let W he a two-parameter Brownian motion on M_|_ x M. In more 
technical terms, W is the orthogonal Gaussian martingale measure on ]R_|_ x R with 
covariance EW{[0, s]x A)W{[0,t]x B) = {s At)Leh{AnB) fors,t G IR+ and bounded 
Borel sets A,B CR. Then 

y{t,r) = aay/K (p^2(^t~s){r - z) dW{s, z) (3.12) 

J J[0,t]xR 

while 

a{t,r) = o-a^/n ip^2s{r - z) dW{s, z). (3.13) 

J J[0,t]xR 

By the equations above we mean equality in distribution of processes. They can 
be verified by a comparison of covariances, as the integrals on the right are also 
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Gaussian processes. Formula (|3.12|) implies that process {y{t,r)} is a weak solution 
of the stochastic heat equation 

yt = yayrr + aaV^W, y{0,r) = (3.14) 

where W is white noise. (See [Ml-) These observations are not used elsewhere in the 
paper. 

Next we record the limits for the quenched mean processes. The four theorems 
that follow require assumption ()2.2j) of stochastic noise and the assumption that the 
annealed probabilities p{0,j) = Emq(0, j) have span 1. This next theorem is the one 
needed for Theorem 12. II for RAP. 

Theorem 3.2. The finite dimensional distributions of processes yn{t,r) converge to 
those ofy{t,r) as n ^ oo. More precisely, for any finite set of points {{tj,rj) : 1 < 
j < k} in M_|_ X M, the vector {xinitj^rf) : I < j < k) converges weakly in M^' to the 
vector {y{tj, rj):l<j<k). 

Observe that property (j3.8|l is easy to understand from the limit. It reflects the 
Markovian property 

E-{x:n = J2 p^{x:^:s = y)E^{xn for s < r, 

y 

and the "homogenization" of the coefficients which converge to Gaussian probabilities 
by the quenched central limit theorem. 

Let us restrict the backward quenched mean process to a single characteristic to 
observe the outcome. This is the source of the first term in the temporal correlations 
()2.19|1 for RAP. The next statement needs no proof, for it is just a particular case of 
the limit in Theorem 13.21 

Corollary 3.3. Fix r G M. As n ^ oo, the finite dimensional distributions of 
the process {yn{t,r) : t > 0} converge to those of the mean zero Gaussian process 
{y{t) : t > 0} with covariance 

Ey{s)yit) = ^{y/TT^ - Vt^) is<t). 
V 27r 

Then the same for the forward processes. 

Theorem 3.4. The finite dimensional distributions of processes an converge to those 
of a as n ^ oo. More precisely, for any finite set of points {{tj,rj) : 1 < j < k} 
in M+ X M, the vector {cLnitj.rj) : I < j < k^ converges weakly in to the vector 
{a{tj,r^):l<j<k). 

When we specialize to a temporal process we also verify path-level tightness and 
hence get weak convergence of the entire process. When r = g in (j2.15|) we get 

Tq{{s At,r),{s At,r)) = CaVs A t 
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with Ca = c"|)/(/5-\/7raf). Since s At is the covariance of standard Brownian motion 
B{-), we get the following limit. 

Corollary 3.5. Fix r G M. As n oo, the process {a„(t,r) : t > 0} converges 
weakly to {B{ca\/t) ■ t > 0} on the path space -Dr[0, oo). 

4. Random walk preliminaries 

In this section we collect some auxiliary results for random walks. The basic 
assumptions, (j2.2|) and span 1 for the p{0,j) = Emo(0, j) walk, are in force throughout 
the remainder of the paper. 

Recall the drift in the ei direction at the origin defined by 

x€Z 

with mean V = —b = K{D). Define the centered drift by 

g{u) = D{uj) -V = - V). 

The variance is ajy = IE[5'^]. The variance of the i.i.d. annealed walk in the Ci direction 
is 

al = Y,i^-V)'Eu-{0,x). 
These variances are connected by 

Let Xn and Xn be two independent walks in a common environment cu, and Yn = 
Xn—Xn- In the annealed sense Yn is a Markov chain on Z with transition probabilities 

g(0, y) = J2 ^K(0, z)uo{0, z + y)\ {y G Z) 
(l{x^y) = ^p{^rz)p{Q,z + y - x) (x 7^ 0,y G Z). 

Yn can be thought of as a symmetric random walk on Z whose transition has been 
perturbed at the origin. The corresponding homogeneous, unperturbed transition 
probabilities are 

q{x, y) = g(0, y - x) = ^p(0, z)p{0, z + y - x) {x,y E Z). 

The g-walk has variance 2cr^ and span 1 as can be deduced from the definition and 
the hypothesis that the p-walk has span 1. Since the g-walk is symmetric, its range 
must be a subgroup of Z. Then span 1 implies that it is irreducible. The g-walk is 
recurrent by the Chung-Fuchs theorem. Elementary arguments extend irreducibility 
and recurrence from q to the g-chain because away from the origin the two walks are 
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the same. Note that assumption (j2.2|) is required here because the g-walk is absorbed 
at the origin iff (j2.2|) fails. 

Note that the functions defined in ()2.7|) are the characteristic functions of these 
transitions: 

X{t) = 5(0' a;)e'*" and \{t) = ^ g(0, x)e''\ 

X X 

Multistep transitions are denoted by q^{x,y) and (f{x,y), defined as usual by 
y) = l{x=y}, q^{x, y) = q{x, y), 

fix,y)= ^ q{x,Xi)q{xi,X2)---qixk-uy) {k>2). 

xi,...,a;fc_ieZ 

Green functions for the q- and g-walks are 

n n 

Gn{.x, y) = Y q^{x, y) and y) = ^ q^{,x, y). 

k=0 k=0 

Gn is symmetric but Gn not necessarily. 

The potential kernel a of the g-walk is defined by 



a{x} = 

n—*oo 



lim{G'„(0,0)-G„(x,0)}. (4.1) 
It satisfies a(0) = 0, the equations 

a{x) = Yq{x,y)a{y) for x 7^ 0, and ^ g(0, ?/)a(y) = 1, (4.2) 



and the limit 



hm ^ = (4.3) 

;^±oo \x\ Zaz 



These facts can be found in Sections 28 and 29 of Spitzer's monograph ^T] . 

Example 4.1. If for some k e Z, p{0,k) + p{0,k + 1) = 1, so that q{0,x) = for 
X ^ {-1, 0, 1}, then a{x) = \x\ /{2al). 

Define the constant 

^ = Y^q{Q^x)-a{x). (4.4) 

x& 

To see that this definition agrees with ()2.8|) . observe that the above equality leads to 

n n 

^ = { E ^'(0' 0) - E E ^(0' ^)^"'(^' 0) }• 

fc=0 fc=0 X 
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Think of the last sum over x as P[Yi + Yk = 0] where Yk is the g-walk, and Yi and 
Yfc are independent. Since Yi + Y^ has characteristic function \(t)\'^{t), we get 

/3 = hm — r (1 - A(t)) y A'^m dt = — r dt. 
Ferrari and Pontes J3] begin their development by showing that 

ViC(s) 

where ( and C are the generating functions 

oo oo 
fc=0 fc=0 

Our development bypasses the generating functions. We begin with the asymp- 
totics of the Green functions. This is the key to all our results, both for RWRE and 
RAP. As already pointed out, without assumption ()2.2|) the result would be com- 
pletely wrong because the g-walk absorbs at 0, while a span h > 1 would appear in 
this limit as an extra factor. 

Lemma 4.1. Let x G M, and let x„ be any sequence of integers such that Xn — n'l'^x 
stays bounded. Then 

hm n-V2G„(x„,0) = ^exp{-f^| dv. (4.5) 

Proof. Por the homogeneous g-walk the local limit theorem ^T] Section 2.5] implies 

that ^ 

hm n-'/''Gn{0,Xn) = ^ I -^exp|-^| dv (4.6) 
n^oo 2ai Jq y/27rv I- 2v J 

and by symmetry the same limit is true for n~^^'^Gn{xn, 0). In particular, 

limn-i/2G„(0,0) = — i=. (4.7) 

n^oo ^T^al 

Next we show 

limn-i/2G„(0,0) = -^. (4.8) 
Using ()4.2j) . a(0) = 0, and q{x,y) = q{x,y) for x 7^ we develop 
g'"(0, x)d{x) = J2 ^"'(O' ^)«(^) = E ^"(0' ^)^(^' 

= J2 ^"(0' ^)^(^' y^^(y) = E ^"^'(O' y^^y^ - ^"(O' o) E ^(O' 
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Identify (3 in the last sum above and sum over m = 0,l,...,r2 — Ito get 
(1 + g(0, 0) + ■ ■ ■ + g"-i(0, 0))/3 = ^ g"(0, x)a(x). 

Write this in the form 

n-i/2^„_^(0, 0)/3 = n-i/^Eo [a(F„)] . 

Recall that Yn = Xn — Xn where Xn and Xn are two independent walks in the 
same environment. Thus by Theorem 13.11 n~^/^y„ converges weakly to a centered 
Gaussian with variance 2(j^. Under the annealed measure the walks Xn and Xn are 
ordinary i.i.d. walks with bounded steps, hence there is enough uniform integrability 
to conclude that vT^^'^Eq \Yn\ — > 2^/ a^jTX. By ()4.Hj) and straightforward estimation, 

1 



n 



This proves (j4.8j) . 

From ()4.7p - ()4.8|) we take the conclusion 



lim ^|/3G'„(0,0) - G„(0,0)| = 0. (4.9) 



Let f{z, 0) = l{^=o} and for A: > 1 let 

f{z, 0) = l{^^o} q{z,Zi)q{zi,Z2)---q{zk-i,0). 

This is the probability that the first visit to the origin occurs at time k, including 
a possible first visit at time 0. Note that this quantity is the same for the q and q 
walks. Now bound 



sup 

zez 



-^Gniz,0) - -^Gn{z,0) 

n \ n 



1 " 

sup ^ V/'^(^,0)|/?G„_fc(0,0) - G„„fc(0,0)|. 



k=0 

To see that the last line vanishes as n — >■ oo, by ()4.9p choose no so that 



|/3G„_fc(0,0) - Gn^kiO,0)\<eV^^ 
for k < n — uq, while trivially 

|/5G„_fc(0,0) - Gn-k{0,0)\<Cno 
for n — riQ < k < n. The conclusion ()4.5p now follows from this and ()4.6|) . □ 
Lemma 4.2. sup sup 0) - + 1, 0)| < oo. 

n>l xeZ 
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Proof. Let Ty = inf{n > 1 : = y} denote the first hitting time of the point y. 

n TyAn n 

Gn{x, 0) = [ ^ l{Yk = 0}] =E,[Y. = 0>] + [ H l^^'^ = 0>" 

fc=0 fc=0 k=TyAn+l 

Ty 

<E,[5];i{n = o}] +G'„(y,o). 

A:=0 

In an irreducible Markov chain the expectation \_ Xlfclo — 0}] is finite for any 
given states x, ?/ jHl Theorem 3 in Section 1.9]. Since this is independent of n, the 
inequalities above show that 

sup sup |G„(x, 0) - ^^(x + 1, 0)1 < oo (4.10) 

n —a<x<a 

for any fixed a. 

Fix a positive integer a larger than the range of the jump kernels q{x, y) and q{x, y). 
Consider x > a. Let a = inf{n > 1 : 1^ < a — 1} and r = inf{?T, > 1 : Yn < a}. Since 
the g- walks starting at x and a; + 1 obey the translation-invariant kernel q until they 
hit the origin, 

PxiYa = y,(T = n]= Px+l[Yr = y + 1,T = Tl]. 

(Any path that starts at x and enters [0, a — 1] at y can be translated by 1 to a 
path that starts at x + 1 and enters [0, a] at y + 1, without changing its probability.) 
Consequently 

n a—1 

Gn{x, 0) - Gn{x + 1,0) = Y^J2P-K = y^^ = ^] {Gn-k{y, 0) - G^-kiv + h 0)) . 

A;=l y=0 

Together with (|4.10|) this shows that the quantity in the statement of the lemma is 
uniformly bounded over x > 0. The same argument works for x < 0. □ 

One can also derive the limit 

lim {G„(0, 0) - G„(x, 0) } = /5~^a(x) 

n— >oo 

but we have no need for this. 

Lastly, a moderate deviation bound for the space-time RWRE with bounded steps. 
Let Xl''^ be the spatial backward walk defined in Section El with the bound ()3.2|) on 
the steps. Let XI' = X*'"^ — i — Vs he the centered walk. 

Lemma 4.3. Form,n G N, let {i{m,n),T{m,n)) G Z^, v{n) > 1, and let s{n) oo 
be a sequence of positive integers. Let a, 7 and c be positive reals. Assume 



'y^^v{n)s{n)" exp{— cg(n)'^} < 00. 



n=l 
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Then for F-almost every u), 



lim max s(n)"P"| max X;^"'")'"^'"'") > cs(n)H7| = q. (4.11) 

n^oo l<m<v{n) ll<k<s(n) J 

Proof. Fix e > 0. By Markov's inequality and translation- invariance, 



P 



uj : max sin) P { max a/ > cs(n)'^^ > > e 

l<m<v{n) Vl<k<s{n) J 



<e^^s{nYv{n)P{ max > cs(n)5+T}. 

l<k<s(n) 

Under the annealed measure P, X^^ is an ordinary homogeneous mean zero ran- 
dom walk with bounded steps. It has a finite moment generating function </'(A) = 
log£'(exp{AX°'°}) that satisfies 0(A) = C(A^) for small A. Apply Doob's inequal- 
ity to the martingale Mk = exp(AX°'° — /c0(A)), note that 0(A) > 0, and choose a 
constant ai such that 0(A) < aiA^ for small A. This gives 

P{ max > cs(n)5+T} < P{ max Mk > exp(cAs(n)5+^ - s(n)0(A))} 

< exp(— cAs(?T.) 2+1' + ais(n)A^) = e"^ ■ exp{— cs(?t.)'^} 

where we took A = s(n)~^. 

The conclusion of the lemma now follows from the hypothesis and Borel-Cantelli. 

□ 



5. Proofs for backward walks in a random environment 

Here are two further notational conventions used in the proofs. The environment 
configuration at a fixed time level is denoted by ujn = {uJx,n : x G Z}. Translations 
on Q are defined by {Tx,n^^)y,k = uJx+y,n+k- 



5.1. Proof of Theorem 13. 2L This proof proceeds in two stages. First in Lemma I3!T] 
convergence is proved for finite-dimensional distributions at a fixed t-level. In the sec- 
ond stage the convergence is extended to multiple t-levels via the natural Markovian 
property that we express in terms of ?/„ next. Abbreviate X^'*'^ = j5^|-"*''J+l-^v^J'l-'^*J_ 
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Then for < s < t, 



z&L 

= E^n^r4-Ln.J = L^^^J + ^}n-'/'{E-{xl:f^^^-^-^) - z) 

z& 

= E^n^S-Ln^J = L^^^J + ^}n-'/'{E-{xl::';^^-^-^) - z) (5.1) 

+ yn{un,r) o Tintb\-[nbu„\,[nt\-[nun\ + n~^/^{[ntb\ - [nsb\ - [nbun\) (5.2) 

where we defined n„ = n~^{\nt\ — \ ns\) so that [nn„J = \nt\ — \ ns\. T^^m denotes 
the translation of the random environment that makes (x, m) the new space-time 
origin, in other words {T^,m^)y,n = ^x+y,m+n- 

The key to making use of the decomposition of yn{t,r) given on fines ()5.ip and 
()5.2|1 is tfiat tfie quencfied expectations 

^ \.^\ns\ I ^'^^ yn\Un,r) O l\ntb\-\nbu.n.\,\nt\-\nur^\ 

are independent because tfiey are functions of environments Um on disjoint sets of 
levels m, wfiile tfie coefficients = \nsh\ + on line ()5.H) converge (in 

probability) to Gaussian probabilities by tfie quencfied CLT as n — > oo. In tfie limit 
tfiis decomposition becomes (|3.8|) . 

Because of tfie little tecfinicality of matcfiing \nt\ — \ns\ witfi \n{t — s)\ we state 
tfie next lemma for a sequence tn ^ t instead of a fixed t. 

Lemma 5.1. Fix t > 0, and finitely many reals ri < r2 < . . . < r^. Let tn he a 
sequence of positive reals such that tn ~^ t. Then as n ^ oo the M.^ -valued vec- 
tor {yn{tn,ri), . . . , ?/„(t„, ttv)) converges weakly to a mean zero Gaussian vector with 
covariance matrix {rq((t,rj), (t,rj)) : 1 < i,j < N} with Tg as defined in ()2.15|) . 

Tfie proof of Lemma I^TTl is tecfinical (martingale CLT and random walk estimates), 
so we postpone it and proceed witfi tfie main development. 

Proof of Theorem \3.^ Tfie argument is inductive on tfie number M of time points 
in tfie finite-dimensional distribution. Tfie induction assumption is tfiat 

[yn{U,rj) ■.l<t<M,l<j<N]^ [y{ti,rj) : 1 < i < M,l < j < N] 

weakly on M*^^ for any M time points < ti < t2 < ■ ■ ■ < and (5-3) 

for any reals ri, . . . , r^r for any finite A^. 
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The case M = 1 comes from Lemma l5.ll To handle the case M + 1, let < ti < 
^2 < ■ ■ ■ < tM+i, and fix an arbitrary (M + l)A^-vector [Oij]. By the Cramer- Wold 
device, it suffices to show the weak convergence of the linear combination 

^ OijUaiU^Tj) = ^ OijyniU^rj) + ^ eM+i,jynitM+i,rj) (5.4) 

l<i<M+l l<i<M i<j<N 

l<j<N i<j<N 

where we separated out the (M + l)-term to be manipulated. The argument will use 
(|5.1|) - (|5.2p to replace the values at tju+i with values at Im plus terms independent of 
the rest. 

For Borel sets B CW define the probability measure 

Apply the decomposition (j5.H) - (j5.2p . with s„ = n~^{\ntM+i\ — ["■^mJ) and 
to get 

z£Z (5.5) 

The 0{n^^/^) term above is n^^^^{yntM+ib\ — \ntMb\ — [?t,s„6J), a deterministic 
quantity. Next we reorganize the sum in (j5.5|) to take advantage of Lemma 15.11 
Given a > 0, define a partition of [—a, a] by 

—a = uq < Ui < ■ ■ ■ < ul = a 

with mesh A = maxju^+i —ui}. For integers z such that —ay/n < z < a^pn, let u{^z) 
denote the value un such that Ujiyfn < z < uij^i^fn. For 1 < j < define an error 
term by 

- {^-(xL;;^^;|J+l«Wv^J'NmJ) _ } j (5.6) 

With this we can rewrite ()5.5|) as 

ynitM+l.Tj) = ^p'^jiUen^^'^, U^+iU^''^ ]yn{tM, Ui) + VniSn, Tj) 



1=0 

^1/4n 



(5.^ 



+ Kj(a) + C(n- 
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Let 7 denote a normal distribution on M with mean zero and variance cr^(tM+i —thi)- 
According to the quenched CLT Theorem 13. H 

p'^ .jiuin^^'^, M^+in^/^ ] 7(m£ — r^, — r^] in P-probabihty as n ^ oo. (5.9) 

In view of fl5.4|l and ()5.8|1 . we can write 

Y dijyniturj) = Pn,i,kyn{ti,Vk) + ^ AI+l,j VniSn, Tj) 

(5.10) 

i?„(a) + 0(n"i/^). 



l<i<M+l l<i<A/ l<i<A' 

l<j<N l<k<K 



Above the spatial points {vk} are a relabeling of {rj, ue}, the cu-dependent coefficients 
Pn,i,k contain constants 9ij, probabilities ,,(M^n^/^, M^+in^^/^ ], and zeroes. The con- 
stant limits Pnik ~^ Pi,k exist in P-probability as n ^ oo. The error in ()5.10|) is 
-Rn(a) = JZj^M+ijRnjio,). 

The variables yn{sn,i^j) are functions of the environments {um '■ [^t-^m+i] > m > 
Mm]} and hence independent of yn{ti,Vk) for 1 < i < M which are functions of 
{um : [ntM] > m > 0}. 

On a probability space on which the limit process {y{t,r)} has been defined, let 
yitu+i — iui ■) be a random function distributed like y{tM+i — tu) •) but independent 

of {yit.r)}. 

Let / be a bounded Lipschitz continuous function on M, with Lipschitz constant 
Cf. The goal is to show that the top line 1)5.111) below vanishes as n ^ oo. Add and 
subtract terms to decompose 1)5.111) into three differences: 

E/( 0:,MU,r,))-Ef(^ Y (^^MU^Tj)) (5.11) 

l<«<Af+l l<i<M+l 
l<j<V i<i<v 



|e/( Y ^^,MU,ri, 

^ l<i<M+l 

i<i<-V 

-E/( Y Pn,,kyniU,Vk)+ Y OM+l,USn,r,))\ (5.12) 



l<i<M l<j<N 
l<k<K 



+ ] ^•^( Y Pn,i,kyn{U,Vk) + Y ^M+l,jyn{Sn,rj)^ 



l<i<M l<j<N 
l<k<K 



Ef(^ Y Pi,kyiti,Vk) + Y ^M+i,jyitM+i-tM,rj)^ I (5.13) 



l<i<M l<j<N 
l<k<K 
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+ ] ^•^( 5Z Pi,ky{ti,Vk) + ^ OM+l,jy{tM+l-tM,rj)^ 
^ l<i<M l<j<N 

l<k<K 

-Efi^ e,,y{U,r,))\. (5.14) 

l<i<A/+l ^ 

The remainder of the proof consists in treating tlie three differences of expectations 

By the Lipschitz assumption and (jS.lOj) . the difference ()5.12|1 is bounded by 

CfE\Rr,{a)\ + 0{n-^/^). 

We need to bound Rn{o). Recall that 7 is an A/'(0, a^itM+i — tA,/))-distribution. 

Lemma 5.2. There exist constants Ci and such that, if a > Oq, then for any 
partition {u^} of [—a, a] with mesh A, and for any 1 < j < N , 

lim sup < Ci(v^ + 7(-cx), -a/2) +7(0/2,00)). 

n— >oo 

We postpone the proof of Lemma 15.21 From this lemma, given 5 > 0, we can 
choose first a large enough and then A small enough so that 

lim sup [ difference (|5.12p ] < e/2. 

n— >oo 

Difference fl5.13|) vanishes as n — > 00, due to the induction assumption ()5.3|) . the 
limits Pnik ^ Pi,k ill probability, and the next lemma. Notice that we are not trying 
the invoke the induction assumption ()5.3|) for M + 1 time points {ti, . . . ,tM,Sn}- 
Instead, the induction assumption is applied to the first sum inside / in fl5.13|) . To 
the second sum apply Lemma I^TTl noting that s.„ ^az+i — tM- The two sums are 
independent of each other, as already observed after ()5.im) . so they converge jointly. 
This point is made precise in the next lemma. 

Lemma 5.3. Fix a positive integer k. For each n, let Vn = {V^, . . . ,V^), Xn = 
(X^, . . . ,X!^), and (n be random variables on a common probability space. Assume 
that Xn and Cn are independent of each other for each n. Let v be a constant k- 
vector, X another random k-vector, and ( a random variable. Assume the weak 
limits Vn V, Xn — > X , and Cn ^ C hold marginally. Then we have the weak limit 

Vn-Xn + Cn ^ V-X + C, 

where the X and ( on the right are independent. 
To prove this lemma, write 



Vn-Xn + Cn= iVn-v)-Xn + V-Xn + C 
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and note that since Vn ^ v in probability, tightness of imphes that (Vn — v) ■ 

Xn — in probabihty. As mentioned, it apphes to show that 

hm [ difference fl5.1Hj) ] = 0. 

n— >oo 

It remains to examine the difference ()5.14j) . From a consideration of how the 
coefficients p^ik (jS-lUj) arise and from the hmit ()5.9p . 



^ pi,kyiU,Vk)+ ^ eM+ijy{tM+i-tM,rj)= ^ Oijy{t 

l<i<M ^<j<N l<i<M 

l<k<K l<j<N 

L-l 

+ ^ (^M+i,j ( X] ^("^ ~ rj,Ui+i - rj]y{tM, ue) + y{tM+i - tn, rj) 
i<j<N e=o 

The first sum after the equahty sign matches all but the [i = M + l)-terms in the 
last sum in ()5.14j) . By virtue of the Markov property in ()3.8|) we can represent the 
variables y{tM+i,rj) in the last sum in ()5.14|) by 



yitu+i.Tj) = / (Pa'i{tM+i-tM){u - rj)y{tM,u) du + y{tM+i - tM,rj). 

Then by the Lipschitz property of / it suffices to show that, for each 1 < j < N, the 
expectation 



^=0 



can be made small by choice of a > and the partition {u^}. This follows from the 
moment bounds on the increments of the y-process and we omit the details. We 
have shown that if a is large enough and then A small enough, 

limsup [ difference fj5.14|) ] <e/2. 

n— >oo 

To summarize, given bounded Lipschitz / and e > 0, by choosing a > large 
enough and the partition {ui\ of [—a, a] fine enough. 



lim sup 



l<i<M+l l<j<A/+l 
l<j<Af l<j<N 



< e. 



This completes the proof of the induction step and thereby the proof of Theorem 
O □ 

It remains to verify the lemmas that were used along the way. 
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Proof of Lemma \5.iA We begin with a calculation. Here it is convenient to use the 
space-time walk X^'"^ = (X^'™, m — k). First observe that 

n-l 



k=0 



n-l 



n-l 



(5.15) 



k=0 

From this, for x, y G Z 

n-l 
k=0 



k=0 



'\{{E-{x:'-) -x}- {i?-(xr) - y})'- 



UJ 



0<k<e<n 



(the cross terms for k < i vanish) 



n-l 



k=0 Z,W(: 

n-l 



E E E P^iK'"" = zjE^iXl'"" = wjP^iX^^'" = ajP^iXl'"' = v} 

k=0 z,w,u,v(^Z^ 

X [g{T^uj)g{TuUj) - g{T^uj)g{TuUj) - g{T;,uj)g{TyUj) + g{T.^u})g{T^,u}) 
(by independence ¥.g{T^uj)g{TuUj) = al,l{z=u}) 



n-l 



k=0 
n-l 

2alY,{Po{Yk = 0} - P._,{n = 0} 



k=0 

= 24(G'„_i(0, 0) - Gn-i{x - y, 0)). 

On the last three lines above, as elsewhere in the paper, we used these conventions: X^ 
and Xk denote walks that are independent in a common environment u, = X^—Xk 
is the difference walk, and Gn{x, y) the Green function of 1^- By Lemma f4. 21 we get 
the inequality 



E 



{{E-iX:'^) -x}- {E-iXy'-) ~y}Y <C\x-y\ 



(5.16) 



valid for all n and all x, ?/ G Z. 
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Turning to i?„j(a) defined in ()5.6|) - (j5.7j) . and utilizing independence, 

\ 1/2 

-\ntM^^z,\ntM\\ 7-1/ -v^L"*a/*J+^.L"*mJn\2 "I \ 



n 



z< — a^pn 
z>ay/n 



' \ntM\ 
z<—a^/n 
z>ay^ 

max \z — \ u{^z)\fn\ |"^^^ 

—a^pn<z<a^pn, 

For the last inequality above we used ()5.16|) . bound ()3.4|) on the variance of the 
quenched mean, and then 

By the choice of u{z), and by the central limit theorem if a > 2\rj\, the limit of 
the bound on E|i?„j(a)| as n ^ oo is C{^/A + 'y{—oo, —a/2) +7(0/2,00)). This 
completes the proof of Lemma 15.21 □ 

Proof of Lemma \5.1\ We drop the subscript from t„ and write simply t. For the main 
part of the proof the only relevant property is that nt„ = 0{n). We point this out 
after the preliminaries. 

We show convergence of the linear combination YliLi^iyn(t,ri) for an arbitrary 
but fixed A^- vector 6 = {61, ... , On)- This in turn will come from a martingale central 
limit theorem. For this proof abbreviate = j5^|-"*''J + l-^»v^J'NJ _ Pq^. 1 < A; < ^nt\ 
define 

TV 

Zn,k = n-'/'Y^'E'^9{Txi_^u;) 

i=l 

SO that by 

[nt] N 

E Zn,k = E ^*^-(^' '^*) + Oin-^^"). 

k=l i=l 

The error is deterministic and comes from the discrepancy ()3.7p in the centering. It 
vanishes in the limit and so can be ignored. 
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A probability of the type P'^{Xl_^ = y) is a function of the environments 

{Oj : [nt\ -k + 2<j< [nt\ } 

while g{Ty^s^) is a function of uog- For a fixed n, {zn,k '■ ^ ^ k < lnt\} are martingale 
differences with respect to the fihration 

l^n,k = (^{^j ■■ [nt\ -k + l<j< [nt\ } il<k< [nt\ ) 

with Unfi equal to the trivial cx-algebra. The goal is to show that Ylk^i ^n,k con- 
verges to a centered Gaussian with variance X]i<i j<Ar ^^^^^((i:, rj), (t, rj)). By the 
Lindeberg-Feller Theorem for martingales, it suffices to check that 

[nt\ 

Y,nzl,\Un,k-i]^ J2 0.9,T,{{t,r,),{t,r,)) (5.17) 

k=l i<i,j<N 

and 

[nt} 

Y,nzlkH\zn,k\ > e}\Un,k-i] (5.18) 

k=l 

in probability, as ra — * oo, for every s > 0. Condition (|5.18p is trivially satisfied 
because \zn,k\ < Cn~^^^ by the boundedness of g. 

The main part of the proof consists of checking ()5.17p . This argument is a gener- 
alization of the proof of Theorem 4.1] where it was done for a nearest-neighbor 
walk. We follow their reasoning for the first part of the proof. Since aj^ = K[g'^] and 
since conditioning on Un,k-i entails integrating out the environments ujint\-k+i, 
one can derive 

[nt] 

fc=l 'i-<i,j<N k=0 

where XI and Xl are two walks independent under the common environment u, 
started at ([ntfej + [rji/nj, \nt\) and ([ntfej + \rj^/n\, \_nt\). 
By (USD 

[nt\-l 

aln'l^ E P{Xl = Xi) r,((t,r,), (t,r,)). (5.19) 

A:=0 

This limit holds if instead of a fixed t on the left we have a sequence t„ t. Conse- 
quently we will have proved ()5.17|) if we show, for each fixed pair (^, j), that 

^-1/2 ^ ^p^^xi = Xi} - P{Xl = Xi}) ^ (5.20) 

k=0 

in P-probability. For the above statement the behavior of t is immaterial as long as 
it stays bounded as n — > 00. 
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Rewrite the expression in (j5.2()|) as 



n 



-1/2 



{P{Xl = Xi I - P{Xl = Xi I W„,o}) 



fc=0 

[nij-l k-l 

A:=0 £=0 
[niJ-1 [ntJ-1 

e=o k=i+i 

\nt\ - 1 

t=o 

where the last hne defines Ri. Check that KRiRm = for i m. Thus it is 
convenient to verify our goal ()5.20|) by checking convergence, in other words by 
showing 



£=0 

= ^ E 
— y 0. 



5^ (p{x^ = Xi I - p{xi, = Xi I 

k=£+l ^ 



2 H 



(5.21) 



For the moment we work on a single term inside the braces in (|5.2H) . for a fixed 
pair k > i. Write Ym = X^^ — for the difference walk. By the Markov property 
of the walks [recall ()2.27|) ] we can write 

P{Xl = Xi\Un,i+,}= J2 P^{Xl = x,Xl = ~x} 

X u'^nti-eix, y - x)u'(^t^_e{x, y - x)P{Yk = | Ye+i = y - y) 

and similarly for the other conditional probability 

P{Xl = Xi\U^,,}= Y P-{Xl = x,Xl=x} 

x,x,y,y&Z 

X E[uf„,j„,(x, y - y - x)]P{Yk = \Ye+i = y - y). 
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Introduce the transition probability q{x, y) of the F-walk. Combine the above de- 
compositions to express the {k, i) term inside the braces in (|5.21|) as 

P{Xl = Xi I Un,i+,} - P{Xl = XI I Ur,e} 
= ^ P-{Xl = x,Xi=x}q'~'~\y-y,0) 

x,x,y,y&Z 

x,x z,w. —M<w<M 

-M<w-z<M 

- E[uf„ij_^(a;, u;)uf„4j_^(x, w - z)]^. (5.22) 

The last sum above uses the finite range M of the jump probabilities. Introduce the 
quantities 

m 

p^(x,x + m)= ^ Mf„ij_^(x,y) = ^ Mf„ij„^(x,?/) 
y.y<m y=—M 

and 

0^(a;, X, z, w;) = p^(x, x + w - z) 

- E [p'^{x, X + w; - 2)] . 

Fix (x, x), consider the sum over z and w on line ()5.22|1 . and continue with a "sum- 
mation by parts" step: 

q''~^~\x -x + z,0) (Mf„tj_£(x, w - 2;) 

z,w: -M<w<M 
-M<w-z<M 

^ - X + ^, 0) (Cr (x, 5, 2;, w;) - C^{x, x,z-l,w- 1)) 

-M<w)<M 
-M<w-z<M 

Yl [q'"~^~\x- x + z,0) - q''-^-\x- x + z + l,0)Y'^{x,x,z,w) 

z,w: -M<w<M 
~M<w-z<M 

2M 

+ -x + z + l, 0)C"(x, X, z, M) 

z=0 



RANDOM AVERAGE PROCESS 35 

-1 

z=-2M-l 

By definition of tlie range M, the last sum above vanislies because Cf{x, x, z, —M — 
1) = 0. Take tliis into consideration, substitute tlie last form above into (j5.22j) and 
sum over k = i + 1, . . . , [nt\ — 1. Define tlie quantity 

^e,n{^) = E (/■"'"'(^,0)-g'=-^-^(a; + l,0)). (5.23) 

k=e+i 

Then the expression in braces in ()5.2ip is represented as 

[nt\-l 

Ri= (^{^^ = I ^-.^+i> - ^{^^ = I 

k=e+i 

= Ep"{X; = x,X; = 5} Y1 Ae,n{x-x + z)C^{x,x,z,w) (5.24) 

x,x z,w. —M<w<M 

-M<w~z<M 

2M lnt\-l 

+ Y,P''{Xl = x,Xj = x}J2 Yl <l''~'~\x - + z + l,0)C^{x,x,z,M) (5.25) 
x,x z=0 k=e+i 

= Re,i + Re,2 

where Ri^i and Re^2 denote the sums on lines (j5.24|) and (j5.25p . 

Recall from ()5.21|) that our goal was to show that Sl=o ^ ^Rj ^ as n — ^ oo. 
We show this separately for Rg i and Ri^2- 

As a function of uj, C^(- ■ ■ ) is a function of Lj\^nt\~£ and hence independent of the 
probabilities on line (|5.24p . Thus we get 

E[i?2 J = Y ^ [P'^iX} = X, XI = £}P"{X; = x', Xi = x'}] 

X Y ^ Ai^nix - X + z)Ai^nix' - X + z') 

-M<w<M -M<w'<M 
-M<w-z<M -M<w'-z'<M 

X E[C'^{x,x,z,w)C^{x',x',z',w')] (5.26) 

Lemma W?2\ implies that y4£ „(x) is uniformly bounded over {£, n, x). Random variable 
('^{x,x, z,w) is mean zero and a function of the environments {uJx,int\-e,^x,int\-i}- 
Consequently the last expectation on line (j5.26|) vanishes unless {x, x} fl {x', x'} ^ 0. 
The sums over z,w,z',w' contribute a constant because of their bounded range. 
Taking all these into consideration, we obtain the bound 

mil] < C{P{XI = XI} + P{Xl = XI} + P{Xl = Xj}). (5.27) 
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By (14. 5p we get the bound 

[nij-l 

J2 HRli] < Cn^^/^ (5.28) 
e=o 

which vanishes as n — oo. 

For the remaining sum i?^^2 observe first that 

C^{x, X, z, M) = M-z)- Enf„,j_,(x, M - z). (5.29) 

Summed over < z < 2M this vanishes, so we can start by rewriting as follows: 



x,x 

2M [ntj-1 

E {q'-''\x-x + z + l,0)-q'-'-\x-x,0))a{x,x,z,M) 

z=0 k=i+l 

2M z 

P"{X^ = X, Xi =x}J2Yl ^^'-(^ -x + m, 0)Q{x, X, z, M) 



z=0 m=0 
2M 



= - E P'^i^i = ^' = ^} E "^^'""^^ -X + m, 0)p^(x, x + M-m) 

x,x m=0 

where we abbreviated on the last line 

p^(x, X + M - m) = p^(x, X + M - m) - Ep^(x, x + M-m). 
Square the last representation for Ri^2i take E-expectation, and note that 
W.[pt{x,x + M - m)pt{x' ,x' + M - m')] =0 

unless X = x'. Thus the reasoning applied to i?£ i can be repeated, and we conclude 

that also n'^ E1=o~^ 0. 

To summarize, we have verified ()5.2H) . thereby ()5.2()|) and condition ()5.17|) for the 
martingale CLT. This completes the proof of Lemma f5. 11 □ 

6. Proofs for forward walks in a random environment 

6.1. Proof of Theorem 13.41 The proof of Theorem 13.41 is organized in the same 
way as the proof of Theorem 13.21 so we restrict ourselves to a few remarks. The 
Markov property reads now (0 < s < t, r G M): 

a„(t,r) = a^{s,r) + ^p-{z[:^J'0 = [r^\ + [nsV\+y} 



X n-V4{i?-(z[;;^|:ij"^^J+^'L"^J) - [rV^\ - y - [nt\V} + n-'/'{[ns\V - [nsV\). 
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This serves as the basis for the inductive proof along time levels, exactly as done in 
the argument following (j5.3|) . 

Lemma 15.11 about the convergence at a fixed t-level applies to an{t, ■) exactly as 
worded. This follows from noting that, up to a trivial difference from integer parts, 
the processes a„(t, ■) and yn(t, ■) are the same. Precisely, if S denotes the P-preserving 
transformation on fl defined by {Slj)x,t = ^^-intb\+x,int\-T, then 

The errors in the inductive argument are treated with the same arguments as used 
in Lemma (5.21 to treat (a). 

6.2. Proof of Corollary 13. 5L We start with a moment bound that will give tightness 
of the processes. 

Lemma 6.1. There exists a constant < C < oo such that, for all n eN, 



Proof. From 

E[{E-g{T,..ou;) - E^giT.o.ouf] = 2al{P[Y^ = 0] - = 0]) 

we get 

P[Y^ = 0] < P[Y^ = 0] foralln>OandxGZ. (6.1) 

Abbreviate Zn = for this proof. E'^{Zn) — nV is a mean- zero martingale with 
increments E'^giTz^ui) relative to the filtration Tin = o'{uk : < k < n}. By the 
Burkholder-Davis-Gundy inequality [7j, 



¥.[{E^{Zn) -nV)*^] < CE 



n-l 



k=0 



Expanding the cube yields four sums 

0<k<n 



■<n 0<ki<k2<n 

C ^[{E-9{Tz,,^))\E-g{T,^u)) 

0<ki<k2<n 

+ C J2 ^[iE^9iTz,,u:))\E-g{T,^u))\E'^g{T,^u)^ 

0<ki<k2<k3<n 
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with a constant C that bounds the number of arrangements of each type. Replacing 
some gf-factors with constant upper bounds simphfies the quantity to this: 



0<k<n 



0<ki<k2<n 



0<fci<fc2<A:3<n 

The expression above is bounded by C(n^/^ + n + n^/^). We show the argument for 
the last sum of triple products. (Same reasoning applies to the first two sums.) It 
utilizes repeatedly independence, Kg{TuUj)g(T^uj) = cr|,l{„=^} for u,v E 1? , and ()6.1|) . 
Fix < /ci < /c2 < ^3 < Let Z'^ denote an independent copy of the walk in the 
same environment lo. 



E 



X = ^}P^{Z'k, = v}E{g{T^u)giT,uj)} 

CE\{E-giT^^u))\E-giT^^u)Yp-{Z,, = ZlJ 
CE {E-g{T2,u))\E-g{T2,u)y 

X Y P^{Zk,+i = u,Z',^^,=v} EP-{Zl^,.^_, = Zl^,^_,} 

(walks Z"^ and Z^ are independent under a common u) 
{E-g{T,^u)Y{E-g{T,^u))' 

X J2 P^{Zk,+i = u,Zl 



< CE 



V} 



PiY', = 0) 



CE 



{E-g{T,^u))\E-g{T,^u)Y]p{Yl_,^^ = 0). 



Now repeat the same step, and ultimately arrive at 

E e[(e-^?(t^,^^))^(e-^,(t^,^^))^(e-^(t^,^..))^ 

0<fci<fe2<fc3<n 

<c Y = Q)PiykU.^i = o)p(nV..-i = 0) 

0<fci<fc2<fc3<r!. 



RANDOM AVERAGE PROCESS 



39 
□ 



By Theorem 8.8 in jl2t Chapter 3], 

E[(a„(t + /i,r) -a„(t,r))'(a„(t,r) -a„(t-/i,r))'] < CY?!'' 

is sufficient for tightness of the processes {a„,(t, r) : t > 0}. The left-hand side above 
is bounded by 

E[(a„(t + h,r)- an{t, r))' ] + E[(a„(t, r) - a„(t - /i, r))' ] . 

Note that if /i < l/{2n) then (a„(t-|-/i, r) —an{t, r)) (an(t, r) —an{t — h, r)) = due to 
the discrete time of the unsealed walks, while ii h > l/{2n) then \n{t + h)\ — \nt\ < 
3nh. Putting these points together shows that tightness will follow from the next 
moment bound. 

Lemma 6.2. There exists a constant < C < oo such that, for all < m < n & N , 
E[({E"(Z°'°) - nV} - {^"(^r) - mV}Y] < C{n - mf/\ 



Proof. The claim reduces to Lemma ffj. II bv restarting the walks at time m. □ 

Convergence of finite-dimensional distributions in Corollarv 13.51 follows from The- 
orem The limiting process a(-) = lim an(-,r) is identified by its covariance 
Ed{s)a(t) = Tg[(s A t,r), {s A t,r)). This completes the proof of Corollarv 13.51 



7. Proofs for the random average process 
This section requires Theorem 13.21 from the space-time RWRE section. 

7.1. Separation of effects. As the form of the limiting process in Theorem 12.11 
suggests, we can separate the fluctuations that come from the initial configuration 
from those created by the dynamics. The quenched means of the RWRE represent 
the latter. We start with the appropriate decomposition. Abbreviate 



Xn,r = x{n,r) = [ny\ + [ry^\. 
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is variable and 



Recall that we are considering y G M fixed, while (t, r) G IR+ x 
serves as the index for the process. 



a 



\,^{xn,r + [nth\) - ao"(x„,.) = ao-(Xp;;' 



X 



x{n,r)-\-\_ntb\ , [ntj 
[ntj 



ntb\ , [ntj 



,x(n,r-)} ^^0(0 
fi=x(n,r)+l 

x(n,r) 

1 f v^x(n,r)+ [ntbj , [ntj , , , "I / 



^x(n,r)+\ntb\, \ nt\ 
%t\ 



j^3:(n,r)+ [n£6J , \ nt\ ^ 



x(n,r)+ [ntbj , \nt\ 



= J2 ^^{^<< 

i>x{n,r) i<x(n,r) 

Recalling the means g{i/n) = E?7q(z) we write this as 

^lnti(^n,r + lntb\) - ao"(x„,,) = Y^{t,r) + H'^{t,r) 

where 

r"(t, r) = Y^ii^ii) - g{r/n)) (l{z > x„,}P-{z < X^i^f^+L-^J' L-J | 



(7.1) 



and 



H''{t,r) = J]f?(z/n)(l{z > x„,jP'^{z < X, 



a;(n,r)+[nt6J , [ntJ ' 
[ntj 



The plan of the proof of Theorem l2.1l is summarized in the next lemma. In the pages 
that follow we then show the finite-dimensional weak convergence n~ 
and the finite-dimensional weak convergence n~^/^y" — > Y for a fixed to. This last 
statement is actually not proved quite in the strength just stated, but the spirit is 
correct. The distributional limit n~^^'^Y"' — *• Y comes from the centered initial incre- 
ments viQ{i) — g{i/n), while a homogenization effect takes place for the coefficients 
P"{^<X^irf^+Ln*''J'L"*J} which converge to limiting deterministic Gaussian probabil- 
ities. Since the initial height functions and the random environments u that drive 
the dynamics are independent, we also get convergence n~^/^(y" + if") ^ Y + H 
with independent terms Y and H. This is exactly the statement of Theorem 12.11 

Lemma 7.1. Let {Qq, J^q, Pq) be a probability space on which are defined independent 
random variables rj and uo with values in some abstract measurable spaces. The mar- 
ginal laws are P for oj and P for rj, and P'^ = ® P is the conditional probability 
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distribution of {oj^rf), given uj. Let H"(co') and Y'^iuj^rf) be M.^ -valued measurable 
functions of{uj,ri). Make assumptions (i)-(ii) below. 

(i) There exists an M.'^ -valued random vector H such that H"(cij) converges weakly 
to H. 

(ii) There exists an M.^ -valued random vector Y such that, for all 9 G M^, 

£"^[6^^-^"] ^ E{e'^-^) m ¥-probabiUty asn^oo. 
Then H" + Y" converges weakly to H + Y, where H and Y are independent. 
Proof. Let 6', A be arbitrary vectors in M^. Then 

< |E [e*"-"" (E'^e^^-^" - Ee''-^)] \ + | (Ee^^'"" - Ee^^'") Ee^^'^l 

< |E [e'^-H" (E'^e^^-Y" - Ee^''^)] \ + |Ee^^-«" - Ee'^-^\ . 

By assumption (i), the second term above goes to 0. By assumption (ii), the integrand 
in the first term goes to in P-probabihty. Therefore by bounded convergence the 
first term goes to as n — oo. □ 

Turning to the work itself, we check first that H^{t,r) can be replaced with a 
quenched RWRE mean. Then the convergence if" H follows from the RWRE 
results. 

Lemma 7.2. For any S,T < oo and for F-almost every uj, 



lim sup n 

n-*oo o<t<T 
-S<r<S 



-1/4 



H-{t, r) - e{y) ■ Vnt\ _ 



0. 



Proof. Decompose H'^(t,r) = H^{t,r) — if^(t,r) where 

H^t, r)= J2 PH^< X^j^f^^^-''^' } . g{^/n), 

i>x(n,r) 

H^{t,r)= j2 pn^><rf^^^"*'^'^"*M-^^(v^)- 

i<x{n,r) 

Working with H^{t,r), we separate out the negligible error. 

H-it,r)= e{y) J2 p-|z<X;;;;)+Lr.*.J'Lr.*j| 



i>x{n,r) 



i>x(n,r) 

Q{y) ■ E' 



f i?i(t,r) 
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rx(n,r)+lntb\, [nt] 



m=l 



— + — ) - ^(y) 

n n / 



Fix a small positive number 5 < |, and use the boundedness of probabilities and the 
function g. 

L„l/2+6j 

mt^r)\< \o{^ + '^)-9{y)\ 



m=l 



c- 



E 



P'^ {Xnr+m< X 



x{n,r)+\ntb\, [nij 
\nt\ 



(7.2) 



m=[ni/2+«J+i 



By the local Holder-continuity of q with exponent 7 > |, the first sum is oiji^f^) if 
(5 > is small enough. Since ^^("'^)+l-"*^J'l-"*J _ ^.^^^ _|_ [^j7,^5J and by time \nt\ the 
walk has displaced by at most M[ntJ, there are at most 0{n) nonzero terms in the 
second sum in (|7.2|) . Consequently this sum is at most 



Cn ■ P"^ {^^in,r)+lntb\, [ntj 



By Lemma 14.31 the last line vanishes uniformly over t G [0,T] and r G [— S", 5] as 
n ^ oo, for P- almost every u. We have shown 



lim sup n-i/^ H'^it, r) - g{y) ■ E'^KXf 



ui\ ( •trx{n,r)+\nth\, \ nt\ 



0<t<T 
-S<r<S 



nt\ 



P-a.s. 



Similarly one shows 



lim sup n 

-S<r<S 



-1/4 



L^l ( yx{n,r)+[ntb}, [nt\ 



nt\ 



■Y] 



P-a.s. 



The conclusion follows from the combination of these two. 



□ 



For a fixed n and the process j^^ i^x^^'^'^^^^^^^'^^^^ — Xn,r) has the same distri- 
bution as the process yn{t,r) defined in (|3.6|) . A combination of Lemma 17.21 and 
Theorem 13 . 21 implv that the finite-dimensional distributions of the processes n~^^^Hn 
converge weakly, as n — * oo, to the finite-dimensional distributions of the mean-zero 
Gaussian process H with covariance 



EH{s,q)Hit,r) = g{yyTg{{s,q),{t,r)). 



(7.3) 
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7.2. Finite-dimensional convergence of Y^. Next we turn to convergence of the 
finite-dimensional distributions of process Y"- in (j7.1|) . Recall that B{t) is standard 
Brownian motion, and cr^ = E[{Xl'^ — V)"^] is the variance of the annealed walk. 
Recall the definition 

POO 

ro((s, g), (t, r)) = / P[aaB{s) >x- q]P[aaB{t) >x-r]dx 

JqVr 

- {l{r>q} I P[(TaB{s) > X - g]P[a„5(t) < X - v] dx 
Jq 

+ J' P[a,B{s) <x- q]P[aaB{t) > x - r] rfx} 

/qAr 
P[aaB{s) <x- q]P[aaB{t) < x - r] dx. 
-oo 

Recall from ()2.9j) that v{y) is the variance of the increments around \ny\. Let 
{Y(t,r) : t > 0, r G M} be a real- valued mean-zero Gaussian process with covariance 

EY{s, q)Y{t, r) = v{y)Vo{{s, q), (t, r)). (7.4) 
Fix and space-time points (ti, ri), . . . , (1^, rjy) G IR+ x R. Define vectors 

Y" = n-i/^ n), . . . , r^)) and Y = (F(ti, n), . . . , Y{t^, r^)) . 

This section is devoted to the proof of the next proposition, after which we finish the 
proof of Theorem 12.11 

Proposition 7.1. For any vector 9 = (6*1, ...,6n) G R^, E'^(e*^''^") ^ E{e^^-^) in 
F-probability as n ^ oo. 

Proof. Let G be a centered Gaussian variable with variance 

N 

S = v{y) ^ Wo((4,rfc),(t;,n)) 

k,l=l 

and so 6* ■ Y is distributed like G. We will show that 

E'^ie'^-^") E{e'^) in P-probability. (7.5) 
Recalling the definition of Y'^{t,r), introduce some notation: 
C(^,t,r) = l{z > x„,.}P-{z < x^(j'-)+Ln*^J'Ln*J| 

- 1{^ < a;„,,}P"{z > x^;;j'')+Lntf>J, NJ I 

so that 

>^'^(t,r) = ^(%"(z)-f,(z/n))C(^,t,r). 



44 MARTON BALAZS, FIRAS RASSOUL-AGHA, AND TIMO SEPPALAINEN 

Then put 

N 
k=l 

and 

U^{i)=n-'''{r,-{^-Q{i/n))u-{z). 

Consequently 

To separate out the relevant terms let 5 > be small and define 

j=[nj/J-[ni/2+iJ 

For fixed lo and n, under the measure P'^ the variables \UrSS)\ are constant mul- 
tiples of centered increments rj^iS) — Q{i/n) and hence independent and mean zero. 
Recall also that second moments of centered increments ?7o (0 ~ Qi^/"^) are uniformly 
bounded. Thus the terms left out of Wn satisfy 

i:\i-[ny\ \ >n^/'^+^ 

and we wish to show that this upper bound vanishes for P- almost every as n oo. 
Using the definition of bounding the sum on the right reduces to bounding 

sums of the two types 

i:\i-[ny\ \ >nl/2+* 

and 

r.\i-lny\ \ >n^/'^+^ 

For large enough n the points a;(n, r^) lie within |?7,^/^+'^ of \ny\ , and then the previous 
sums are bounded by the sums 

^-1/2 ^ ^pa;|^ ^ j5^|r(n,r-fc)+LnffcbJ, Ln*feJ j 

i>a;(n,rfc)+(l/2)nl/2+« 

and 

^-1/2 ^ ^pi^l^ ^ jl^a;(n,rfc)+[ntfc6J, [ntfej I 

i<x(n,rfc)-(l/2)nl/2+i 
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These vanish for P-almost every as n oo by Lemma 14 .Hf in a manner similar 
to the second sum in erg). Thus 'E'^[{Wn - ■ Y")2] and our goal has 
simplified to 

E'^(e*^") ^ E{e'^) in P-probability. (7.6) 

We use the Lindeberg-Feller theorem to formulate conditions for a central limit 
theorem for Wn under a fixed u. For Lindeberg-Feller we need to check two conditions: 

[nyJ + [ni/2+'5j 

(LF-i) ^ V [[/„(z)2] S 

j=[n5j-[ni/2+«J 

(LF-ii) ^ E'^ [f/„(i)2 . l{lu„i^)\>e}] for all e > 0. 

i=[nyJ-[nl/2+iJ 

To see that (LF-ii) holds, pick conjugate exponents p,q > 1 {l/p+ l/q = 1): 

E- [Unit)' ■ liU.(^)^>e^}] < (E^ (P^ [Un{tf > s'])' 

< £"f (E- [lUrM'"])'' (E- [Unit)'])' < Cn-^/2-i/{2,)_ 

In the last step we used the bound |t/n(^)| < Cn~^^^ IVoi'^) ~ Q{'i'/^)\^ boundedness of 
Q, and we took p close enough to 1 to apply assumption (|2.1ip . Condition (LF-ii) 
follows if (5 < l/(2g). 

We turn to condition (LF-i). 

S-{u)= ^^[Unm= E n-"M^ln)[v-n{i)f 

i=[nyJ-[ni/2+«J i=[n5j-[nl/2+iJ 
Ln5j + [ni/2+«J LnyJ + [nl/2+^J 

5^ n"V2b(z/n)-i;(y)]K(z)]^+ 5^ r."i/2t;(y)K(z)]2 

j=[n5j-[nV2+iJ j=|^„j^j_|^„i/2+«j 

Due to the local Holder-property (j2.1(jp of the first sum on the last line is bounded 
above by 

for sufficiently small 5. Denote the remaining relevant part by S^{u)), given by 

LnyJ + Lni/2+«J [nV2+ij 

5"(^)= 5^ n"i/2.;(y)K(z)]2 = .;(y)n-i/^ ^ «(m+HJ))' 

j=[n5j-[ni/2+«J m=-Lnl/2+iJ 
AT [ni/2+«J 

= z;(y) ^ ^ Cn{[ny\+m,tk,rj,)Q{[ny\+m,tuTi). (7.7) 
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Consider for the moment a particular [k, I) term in the first sum on hne (j7.7|) . Rename 
(s, g) = {ti:,rk) and (t,r) = {ti,ri). Expanding the product of the C^'^^ctors gives 
three sums: 

L„l/2+5j 

m=-[ni/2+^J 
L„l/2+ij 

m=-[ni/2+<5j 

X p-(x^^J'-)+L-*^J'L-*J < LnyJ +m) 

-Li 1 pw/ T^a:{n,g)+[ns6J, [nsj , , , \ 

+ -'-{m,<Lgv^J}-'-{m>Lrv/HJ}-^ [^lns\ < [n?/J + mj 

X ^X'M+Intbi, lnt\ ^ ^^-j ^ I ^^_g^ 

L„l/2+ij 

+ n ' l{m<L5x/HJ}l{m<Lrv^J}P 1,-^ < [nt/J + mj 

m=-[ni/2+«J 

Each of these three sums ()7.8p - (j7.10p converges to a corresponding integral in P- 
probability, due to the quenched CLT Theorem 13.11 To see the correct limit, just 
note that 

and recall that — 6 = \^ is the average speed of the walks. We give technical details 
of the argument for the first sum in the next lemma. 

Lemma 7.3. As n oo, the sum in ()7.8|) converges in F-probability to 

POO 

/ P[aaB{s) > X - q]P[aaB{t) > X -r]dx. 

J q\Jr 

Proof of Lemma \7J^ With 
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and 



f:{x)dx, 



qyr 



the sum in ()7.8|) equals l':^ + 0{n ^/^). By the quenched invariance principle Theorem 
13. H for any fixed x, fn{x) converges in P-probability to 

fix) = P[aaB{s) >x- q]P[aaB{t) > x - r]. 

We cannot claim this convergence P- almost surely because the walks ^|^^^y'')+l-"*''J' L'^^J 
change as n changes. But by a textbook characterization of convergence in prob- 
ability, for a fixed x each subsequence n(j) has a further subsequence n{i() such 
that 



P 



^ • fn{j,){x) fix) 



1. 



By the diagonal trick, one can find one subsequence for all x G Q and thus 



VHj)},3{j,}:P ^ : Vx G Q : 



Of) 



x 



fix) 



Since and / are nonnegative and nonincreasing, and / is continuous and decreases 
to 0, the convergence works for all x and is uniform on [q V r, oo). That is. 



V{n(j)},3{ja:P 



UJ : 



I f 



f\ 



L°°[(jVr, oo) 







It remains to make the step to the convergence of the integral to j^^, f{x) dx. 
Define now 



J^iA) 



fnix)dx. 



qVr 



Then, for any A < oo 



f{x)dx 



q\/r 



1. 



In other words, JniA) converges to J^^^ f{x)dx in P-probability. Thus, for each 
< A < oo, there is an integer m{A) such that for all n > m{A) 



P 



■KiA) 



f{x)dx 



qyr 



>A-^ 



<A-\ 



Pick An oo such that m{An) < n. Under the annealed measure P, X'^^ is a 
homogeneous mean zero random walk with variance 0{n). Consequently 



E[|/--J-(A„,)|]< 



A„An* 



nf:ix)]dx 



< 



jji ^x{n,r) + \ntb\,\nt\ 



> x{n, r) — \ry/n\ + \_x\/n\ j dx — > 
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Combine this with 



P 



LU : 



f{x)dx 



<A-\ 



Since J^"^ f{x)dx converges to f{x)dx, we have shown that converges to this 
same integral in P-probabihty. This completes the proof of Lemma 17. HI □ 

We return to the main development, the proof of Proposition 17.11 Apply the 
argument of the lemma to the three sums ()7.8j) - ()7.1()j) to conclude the following limit 
in P-probability. 

lim n"^/^ V Cnilny\ +m,s,q)Cn{lny\ +m,t,r) 

n^oo ' ' 

m=-[ni/2+«J 

/■oo 

= / P[(Ta5(s) > X - q\P\aaB{t) > x-r\dx 

J q\jT 

- |l|,>,} I P[(JaB{s) > X - q]P[aaB{t) <x-r]dx 

Jq 

+ l{q>r} P[(raB{s) < X - q]P[aaB{t) >x-r] c/x} 

/gAr 
PKB{s) < X - q]P[aaB{t) <x-r]dx 
■oo 

= To{{s,q),{t,r)). 

Return to condition (LF-i) of the Lindeberg-Feller theorem and the definition ()7.7|) 
of S"'{uj). Since S"'{uj) — S"'{u) — >• as pointed out above (|7.7|1 . we have shown that 
S*" — > S* in P-probability. Consequently 



1. 



This can be rephrased as: given any subsequence {n{j)}, there exists a further sub- 
sequence {n{je)} along which conditions (LF-i) and (LF-ii) of the Lindeberg-Feller 
theorem are satisfied for the array 

{Unud^) ■■ [nij,)y\ - [n{jeY/''-'\ < i < yn{3,)y\ + \n{3,fl'^'\ , £ > 1} 

under the measure P"^ for P-a.e. uo. This implies that 



V{n(j)},3{ja:P 



E{e 



1. 



But the last statement characterizes convergence E'^(e^^") ^ E{e'^) in P-probability. 
As we already showed above that Wn — ■ Y" — > in P'^-probability P-almost surely, 
this completes the proof of Proposition 17.11 □ 



RANDOM AVERAGE PROCESS 



49 



7.3. Proofs of Theorem 12.11 and Proposition 12.21 

Proof of Theorem \2. 1\ The decomposition ()7.1|) gives = n^^^^(Y^ + H^). The 
paragraph that follows Lemma 17.21 and Proposition 17.11 verify the hypotheses of 
Lemma 17.11 for iJ" and F". Thus we have the limit z = Y + H in the 

sense of convergence of finite-dimensional distributions. Since Y and H are mutually 
independent mean-zero Gaussian processes, their covariances in fl7.3|l and ()7.4|1 can 
be added to give (jTTT^ . □ 

Proof of Proposition \2.'A The value ()2.23|) for (3 can be computed from ()2.8|) . or from 
the probabilistic characterization ()4.4j) of (3 via Example 14.11 If we let u denote a 
random variable distributed like Uo(0, —1), then we get 

Em - E(m2) _ E(m2) - (Em)2 

^ ~ Em - (Em) 2 " ~ Em - E(m2) ' 

With obvious notational simplifications, the evolution step ()2.22|) rewrites as 

7]'{k) - p = (1 - Uk){r]{k) - p) + Uk^i{r]{k - 1) - p) + (m^^i - Uk)p. 

Square both sides, take expectations, use the independence of all variables {'ri{k — 
1), ?7(/c), Mfc, Mfc_i} on the right, and use the requirement that rj'{k) have the same 
variance v as rj{k) and ri{k — 1). The result is the identity 

v = v{l- 2Em + 2E(m2)) + 2p2(E(M^) - (Em)^) 

from which follows v = hip^. The rest of part (b) is a straightforward specialization 
of (jnUl- □ 
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